MORE ACCURATE APPROXIMATIONS 
FOR THE GAMMA FUNCTION 
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Abstract. A series transformation idea inspired by a formula of R. W. Gosper 
and some asymptotic expansions for the central binomial coefficients leads us 
to new accurate approximations for the Gamma function. 



1. Introduction 

The Gamma function plays an important role in several fields of mathematics 
such as probability theory or combinatorics. One often has to evaluate the function 
for large positive values. One way to aim this is to use asymptotic approximations. 
It is well known that for large values of x the Gamma function has the asymptotic 
series of the form [3l O HI [7] 
(1.1) 

1 1 139 571 



r(x + i) 



I2x 288x2 51840x3 2488320x4 



Equation (jl.ip is called Stirling's formula however, Laplace was the first who derived 
it by his approximation method for special integrals. Another famous result is the 
Stirling series [31 El [7] 

(1.2) logr(x + l)^ (x + i)log^-^ + ilog(27r) + -^^ ^ ^ 



2j ^ 2 ^ ' ' 12x 360x3 1260x5 

A main advantage of this latter series is that it has only odd powers of the variable. 
For the past almost three hundred years several authors established fascinating new 
asymptotic formulas to improve the accuracy of (jl.ip . For example, Karatsuba 
showed that a formula of Ramanujan can turn into an asymptotic expansion 



(1.3) r (x + 1) ^ x-e--^A^8x3 + 4x2 + x + ^ - ^ + • • •. 

Mortici proved in his more recent paper |lj the following expansion similar to Karat- 
suba's: 



11 31 1 3Q 

(1.4) 7l! = ^(n + l)-n"e-"^AJ2r^+- + 



3 36n 3240^2 77760^3 

We develop some new variants of (|l.ip in this paper and show that these formulas 
are numerically more efficient than much of the early ones in many cases. The 
first few values of the newly introduced coefficients and sequences can be found in 
Appendix El 
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2. New asymptotic expansions 
The motivating examples are the fohowing two asymptotic expansions [7]: 

, , f2n\ 2^" / 1 1 5 21 

(2-1) L l-7- + T^ + 



n 



^ 8n 128^2 10247i3 32768^4 
f2n\ 2^" / 1 21 



V"/ {n + l)\ 64(n+i)' 8192(n+i)' 



(2.2) 



The first one is the standard asymptotic series of the central binomial coefficients. 
If one expands them into a series in powers of 1/ (n + 1/4), the asymptotic series 
contains only even powers. This remarkable result suggests that there might have 
been an asymptotic expansion similar to (|l.ip that is, it contains only even powers 
of the shifted variable. The formula 

(2.3) r(x + i) = 

known as Gosper's approximation 6 can be a good starting point. Our aim is to 
elaborate the asymptotic series part in Gosper's formula. It seems from (|2.2[) that 
another series in terms of 1/ (x + 1/6) would be the right choice. It can be shown 
that a series like that contains even and odd powers. If we insist to have even 
powers only we are lead to the form 
(2.4) 




r(a; + l) -x'^e 




92 



[x + Vif {x + V2)'^ 



where the sequences {ffn}„>o ^^nd {fn}„>o to be determined. One of our main 
result is 

Theorem 2.1. The Gamma Junction has an asymptotic series expansion of the 
form 

(2.5) T{x + l)^x='e- 




/ ^ / I \2n' 



as X oo, where the sequences {5n}„>o '"^'^ {^"}n>o ^^''^ he found from the recur- 
rence 

Here the an coefficients are those appearing in (jl.ll) . i.e., 

(2.7) ^(^+lL^y^. 

n>0 

Proof. As X — )■ 00 we have 



T{x + l)e- _ r(x+l)e- / ^ 1 \-^^\ + 
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From the binomial formula wc find 



1 



1 

6a; 



-1/2 



E 



-1/2 



n 1 6"a;" 



as a; — ?> cx). Thus we obtain the asymptotic expansion 
r (x + 1) e 



Suppose the expansion of the form 



l/2\ o„ 

i 



6^ 




9n 



„>o + 

Now we expand the right hand side in powers of Xjx: 

'-2n\ vL 



E 

n>0 



9n 



(x + Vn) 



E 9n f Vn 
™2n V T 



n>0 



EE 



-2n 



9nVn 




n~2j 



n>Q 1>Q 

According to the uniqueness of asymptotic series, the proof is complete. 
From (ESI) we find 



□ 



1 _ 3857 _ 20932906335329 
go - 1, 51 - — , .92 - " 3210400' ~ 34283052002304000' 



23 



V2 



1792627 



«3 



570984637359867601981 



90' 7289730 " 2288928529497568067550' 

Note that vq vanishes in (|2.5I) due to the zero power. We have obtained an expansion 
in even powers however, the shift sequence {wn}„>g has different terms whereas 
has constant (= 1/4) shift in all terms. Numerical evaluation of the first few v„ 

vi = 0.2555555555 . 
V2 0.2459113026. 
V3 = 0.2494549873 . 
V4 = 0.2498398924 . 
1^5 = 0.2499584970. 
V6 = 0.2499884146 . 
V7 = 0.2499963289 . 



leads us to the 
Conjecture 2.1. lim 



1/4. 



This conjecture suggests a new asymptotic series to the Gamma function in 
terms of 1/ (a; + 1/4). Our second result is 
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Theorem 2.2. The Gamma function has an asymptotic series expansion of the 
form 



(2.9) r (a: + 1) ^ x^e'^ 2n U+^^Y. T^V' 

V ^ ^ fc>o + a) 



as X ^ oo, where the the Gk coefficients are given by 
The a„ coefficients again, are from (fLT|) . 

Proof. Similary to the previous proof we expand our series in terms of 1/x: 
Gk ^ Gfe ^ f-k\ 1 1 



E^k _ 



k>0 + i) fc>0 i>0 



X'" fr-f \ / / 4' a;' 




According to (|2.8p and the uniqueness of asymptotic series the proof of (|2.10p is 
complete. □ 

3. Numerical comparisons 

We win compare in this paragraph the numerical performance of some asymp- 
totic formulas to the Gamma function with our new formulas for large values. We 
compare the following approximation formulas for T {x + 1). 



(Stirling) , 





1 


1 


12^ " 


360x3 


^ 1260x5 


1 


1 


139 




288x2 


51840x3 


1 




1 


^ 2^ " 




240x3 



/x\^ / 1 1 139 \ 

(3.2) ( - V27rx 1 H h ^ ^ (Laplace) , 

^ ' \e) \ 12x 288x2 51840x3 \ v >^ 



(3.3) (£)'v&-/ l + + .- (IU„,a™ja„) 



(3.4) + + (Mortid) 



(3.5) f-V J27r(x+-) 1 + ^ ^ (New). 

^ ' ^e) \j \ q)\ 144(x+i)' 12960 (x+i)' ) 

The following table displays the number of exact decimal digits (edd) of the 
formulas for some values of x. Exact decimal digits are defined as follows: 

approximation (x) 



(3-6) edd (x) = - logio 1 - 

1 (x + 1) 

In the table below the (i)-th entry (i = 1,2, . . .) in a line starting with "name" 
is the edd of the given approximation using the series up to the z-th order term. 
The "— " sign indicates that the approximation is smaller and the "+" sign (not 
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displayed) indicates that the approximation is larger than the true value. Note that 
in the case of Stirling's formula the first n terms of the asymptotic series give the 
2nth order approximation. 



Formula 


X 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Stirling 


100 




8.6 




-13.1 




17.2 




-21.1 


Laplace 


100 


-6.5 


8.6 


11.7 


-13.1 


16.2 


17.2 


-20.4 


-21.1 


Ramanujan 


100 


-5.7 


-9.2 


11.0 


-13.3 


-15.4 


17.3 


19.5 


-21.1 


Mortici 


100 


-6.2 


8.6 


11.4 


-13.1 


-15.9 


17.2 


-20.0 


-21.1 


New 


100 


-6.2 


10.1 


10.9 


-14.9 


-15.2 


19.4 


19.2 


-23.0 


Stirling 


1000 




11.6 




-18.1 




24.2 




-30.1 


Laplace 


1000 


-8.5 


11.6 


15.6 


-18.1 


22.2 


24.2 


-28.3 


-30.1 


Ramanujan 


1000 


-7.7 


-12.2 


15.0 


-18.3 


-21.4 


24.3 


27.5 


-30.1 


Mortici 


1000 


-8.2 


11.6 


15.4 


-18.1 


-21.9 


24.2 


-28.0 


-30.1 


New 


1000 


-8.2 


13.1 


14.9 


-19.7 


-21.2 


26.9 


27.2 


-33.5 


Stirling 


10000 




14.6 




-23.1 




31.2 




-39.1 


Laplace 


10000 


-10.5 


14.6 


19.6 


-23.1 


28.2 


31.2 


-36.3 


-39.1 


Ramanujan 


10000 


-9.7 


-15.2 


19.0 


-23.3 


-27.4 


31.3 


35.5 


-39.1 


Mortici 


10000 


-10.2 


14.6 


19.3 


-23.1 


-27.9 


31.2 


-36.0 


-39.1 


New 


10000 


-10.2 


16.1 


18.9 


-24.7 


-27.2 


33.7 


35.2 


-42.1 



Table 1 . The number of exact decimal digits of the asymptotic 
series for some values of x. 



Conclusion. It is seen that when we use odd order approximations, Laplace's 
formula is the most accurate. In the case of even orders Ramanujan's approximation 
is better than Stirling's, Laplace's and the one by Mortici, but our new formula 
gives better approximations even than that of Ramanujan's. 

Expression (j2.5l) is a slightly different than the previous ones, thus we consider 
it's edds in a separate table. The notations are the same except the fact that this 
series contains only even order terms. 



Formula x (2) (4) (6) (8) (10) 



Special 100 10.9 -15.2 19.2 -22.9 26.5 
Special 1000 14.9 -21.2 27.2 -32.9 38.5 
Special 10000 18.9 -27.2 35.2 -42.9 -50.5 

Table 2. The number of exact decimal digits of the special as- 
ymptotic series ()2.5p for some values of x. 
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Appendix A. Tables of coefficients 

The following table gives the first few values of the sequences {5n}„>o and 
{■^«}n>o appearing in ([23]) . 



50 


1 


vo 




91 


1 

144 


Vl 


23 
90 


92 


3857 


V2 


1792627 


3110400 


7289730 


93 


20932906335329 


V3 


570984637359867601981 


34283052002304000 


2288928529497568067550 



Table 3. 

The table 
and {v„}„>o 



50 


1.000000000000000 


vq 




51 


0.006944444444444 


Vl 


0.255555555555555 


52 


-0.001240033436214 


V2 


0.245911302613402 


53 


0.000610590513759 


V3 


0.249454987345193 


54 


-0.000655407405149 


va 


0.249839892410196 


55 


0.001199164540953 


V5 


0.249958497082160 



Table 4. The first few numerical values of the sequences {5ra}„>o 
and {w„}„>o- 

The following table gives the first few coefficients appearing in the asymptotic 
series (|2J|) . 



Go 


1 


Gs 


53 


Gio 


360182239526821 


2612736 


300361133850624000 


Gi 





Ge 


5741173 


Gil 


104939254406053 


9405849600 


210853515963138048000 


Gs 


1 


G7 


37529 


Gl2 


508096766056991140541 


144 


18811699200 


151814531493459394560000 


G3 


1 


Gs 


710165119 


Gi3 


70637580369737593 


12960 


1083553873920 


151814531493459394560000 


G4 


257 


Go 


3376971533 


Gi4 


289375690552473442964467 


207360 


4022693756928000 


21861292535058152816640000 



Table 5. The first few Gk coefficients. 
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